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Abstract. The pathway model of Mathai (2005) is shown to be inferable 
from the maximization of a certain generalized entropy measure. This en- 
tropy is a variant of the generalized entropy of order a, considered in Mathai 
and Rathie (1975), and it is also associated with Shannon, Boltzmann-Gibbs, 
Renyi, Tsallis, and Havrda-Charvat entropies. The generalized entropy mea- 
sure introduced here is also shown to have interesting statistical properties 
and it can be given probabilistic interpretations in terms of inaccuracy mea- 
sure, expected value, and information content in a scheme. Particular cases 
of the pathway model are shown to be Tsallis statistics (Tsallis, 1988) and 
superstatistics introduced by Beck and Cohen (2003). The pathway model's 
connection to fractional calculus is illustrated by considering a fractional 
reaction equation. 



1 Introduction 

The fundamental problem pursued in equilibrium statistical mechanics is 
that given a large number of physical species, such as atoms, one wishes to 
know how they distribute according to some common property, e.g. veloc- 
ity or energy (Ebeling and Sokolov, 2005). A simple mathematical model 
to understand the problem is at the center of statistics and probability the- 
ory. In order to deal with applications to physical situations of interest 
one takes into consideration the fundamental hypothesis of equal a priori 
probabilities for regions in phase space of an isolated system. This hypoth- 
esis is based on our insufficient knowledge for a specification of the precise 
state of the physical system under consideration. This hypothesis allows us 
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to assign systems to states that agree equally well with our knowledge of 
the actual condition of the system. This leads to the Boltzmann-Gibbs en- 
tropy, or Boltzmann principle as Einstein called it, S — k InW, where W 
is the thermodynamic probability which is defined as the total number of 
equally probable microstates corresponding to the given macrostate. The 
Boltzmann constant is denoted by k. The Boltzmann-Gibbs entropy is rele- 
vant for situations such that all possible states of the system are considered 
equally probable. If we consider such a system in contact with a thermostat 
then we obtain the usual Maxwell-Boltzmann distribution for the possible 
states by maximizing the Boltzmann-Gibbs entropy S with the normaliza- 
tion and energy constraints. However, in nature many systems show distri- 
butions which differ from the Maxwell-Boltzmann distribution. These are 
usually systems with strong autocorrelations preventing the convergence to 
the Maxwell-Boltzmann distribution in the sense of the central-limit theo- 
rem. Well known examples in physics are: self gravitating systems, charged 
plasmas, Brownian particles in the presence of driving forces, and, more gen- 
erally, non-equilibrium states of physical systems (Abe and Okamoto, 2001; 
Gell-Mann and Tsallis, 2004). Then it is natural to ask the question of 
whether non-Maxwell-Boltzmannian distributions can also be obtained from 
a corresponding maximum entropy principle, considering a generalized form 
for the entropy. For this purpose, different forms were proposed, as for in- 
stance the Tsallis entropy S q = w x ^_~ x , where q is the entropic index, that 
is considered the basis for a generalization of Boltzmann-Gibbs statistical 
mechanics (Abe and Okamoto, 2001; Gell-Mann and Tsallis, 2004). In the 
present paper we are investigating the link between entropic functionals and 
the corresponding families of distributions in Mathai's pathway model. We 
come to the conclusion that this link is also important to physically analyze 
fractional reaction equations in terms of probability theory. 

The structure of the paper is the following: In Section 2 we introduce 
basic notions of Mathai's pathway model in terms of parametric families of 
distributions. In Section 3 we introduce a generalized entropic measure and 
investigate its characteristics and establish the link to parametric families of 
distributions in Mathai's pathway model, including Tsallis' distribution. In 
Section 4 we establish the link between a fractional reaction equation, its 
reaction coefficient considered a random variable, and Tsallis statistics and 
superstatistics. 



2 



2 Preliminaries for Mathai's pathway model 



For practical purposes of analysing data of physical experiments and in build- 
ing up models in statistical physics, we frequently select a member from a 
parametric family of distributions. It is often found that fitting experimental 
data needs a model with a thicker or thinner tail than the ones available from 
the parametric family, or a situation of right tail cut off (Honerkamp, 1994). 
The experimental data reveal that the underlying distribution is in between 
two parametric families of distributions. This observation either apeals to 
the form of the entropic functional or to the representation by a distribution 
function. In order to create a pathway from one functional form to another a 
pathway parameter is introduced and a pathway model is created in Mathai 
(2005). This model enables one to proceed from a generalized type-1 beta 
model to a generalized type-2 beta model to a generalized gamma model 
when the variable is restricted to be positive. More families are available 
when the variable is allowed to vary over the real line. Mathai (2005) deals 
mainly with rectangular matrix-variate distributions and the scalar case is 
a particular case there. For the real scalar case the pathway model is the 
following: 

f{x) = cx^~ l [l~a(l-a)x s ]^, (1) 

a > 0, 5 > 0, 1 — a(l — a)x s > 0, 7 > where c is the normalizing constant and 
a is the pathway parameter. For a < 1 the model remains as a generalized 
type-1 beta model in the real case. For a — 1, 7 = 1, 5 — 1 we have Tsallis 
statistics for a < 1 (Tsallis, 1988, 2004). Other cases available are the regular 
type-1 beta density, Pareto density, power function, triangular and related 
models. Observe that (1) is a model with the right tail cut off. When a > 1 
we may write 1 — a = — (a — l),a > 1 so that f(x) assumes the form, 

f(x) =cx^- 1 [l + a(a-l)x s ]-^, x>0 (2) 

which is a generalized type-2 beta model for real x. Beck and Cohen's super- 
statistics belong to this case (2) (Beck and Cohen, 2003; Beck, 2006). For 
7 = l,a = 1,(5 = 1 we have Tsallis statistics for a > 1 from (2). Other 
standard distributions coming from this model are the regular type-2 beta, 
the F-distribution, Levi models and related models. When a — > 1 the forms 
in (1) and (2) reduce to 

f(x) = cx"- x e~ ax \ x > 0. (3) 
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This includes generalized gamma, gamma, exponential, chisquare, Weibull, 
Maxwell-Boltzmann, Rayleigh, and related models (Mathai, 1993a; Hon- 
erkamp, 1994). If x is replaced by \x\ in (1) then more families of distri- 
butions are covered in (1). The normalizing constant c for the three cases 
are available by putting u = a(l — a)x s for a < 1, u — a(a — l)x s for 
a > 1, u — ax 6 for a — > 1 and then integrating with the help of a type- 
1 beta integral, type-2 beta integral and gamma integral respectively. The 
value of c is the following: 

#-«)] ] r(l + ^ + i) . 

for a < 1 



rfflr^ + i 

S[a(a - l)]iT 



1 

, for > 0, a > 1 



r 



, for a — > 1. (4) 



<5 a« 



Observe that in (2) and (3), - also belongs to the same family of densities 
and hence in (2) and (3) one could have also taken x~ 5 with 5 > 0. 



3 Pathway model from a generalized entropy 
measure 

We introduce a generalized entropy measure here. This is a generalization of 
Shannon entropy and it is also a variant of the generalized entropy of order 
a in Mathai and Rathie (1975, 1976). Let us take the discrete case first. 
Consider a multinomial population P = ...,p k ),pi > 0, i — 1, . . . , k, pi + 
. . . +Pk — 1- Define the function 

M Ka {P) = E ^ =lR2 a^l, -oo<a<2 (5) 

a — 1 

k 

limM M (P) = -^m^^P) (6) 

i=i 

by using L'Hospital's rule. In this notation OlnO is taken as zero when any 
Pi = 0. Thus (5) is a generalization of Shannon entropy Sk(P) as seen from 
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(6). Note that (5) is a variant of Havrda-Charvat entropy Hk, a {P) and Tsallis 
entropy Tfc )Q (P) where 

HkAP) = E ltf_\ 1 i a^l, a>0 (7) 

and 

r*,a(P) = E ^ lP? ~ 1 , a^l, a>0. (8) 
1 — a 

We will introduce another measure associated with (5) and parallel to Renyi 
entropy Rk, a in the following form: 

In p 2 ~ a ) 

M* k jP)= V ,=1 1 y , a^l,-oo<a<2. (9) 
a — 1 

Renyi entropy is given by 



In (Ei-iP?) 



a>0. (10) 

1 — a 

It will be seen later that the form in (5) is amenable to power law, pathway 
model etc. First we look into some basic properties enjoyed by Mfc jCt (P). 



3.1 Properties 

(i) No n- negativity: M k a (P) > with equality only when one Pi = 1 and 
the rest zeros. 

(ii) Expansibility or zero-indifferent: Mfc +l Q ,(P, 0) = M^ a (P). If an im- 
possible event is incorporated into the scheme, that is, Pk+i = it will 
not change the value of the entropy measure. 

(iii) Symmetry: Mk, a {P) is a symmetric function of pi,...,Pk- Arbitrary 
permutations of p±, . . . ,pk will not alter the value of Mk, a (P). 

(iv) Continuity: M^ a {P) is a continuous function of pi > 0, i — 1, . . . , k. 

(v) Monotonicity: M kjQ , ( |, . . . , |J is a monotonic increasing function of k. 



5 



(vi) Inequality: M M (pi, ...,p k )< M Ka (|, . . . , |) . 

(vii) Branching principle or recursivity: 

M fc)Q (pi,...,p fe ) = M fe _i )Q (pi + p 2 ,P3, • • • ,Pk) + (pi +P2) 2 ~ a x 
x M 2 , f Pl K ^ 



P1+P2 Pi +P2/ ' 

This property indicates what happens to the measure if two of the mu- 
tually exclusive and totally exhaustive events, defining the multinomial 
population, are combined. 

(viii) Non-additivity: Consider independent multinomial populations P = 
(pi, . . . ,p„) and Q = (<?i, . . . , q m ) such that £™ =1 £Jli p^ = 1, £™ =1 Pi = 
1, Y^jLiQj — 1- Then the joint density is of the form (pi^i,..., 
Piq m , ■ ■ ■ iPnQii ■ ■ ■ ,PnQm)- Let us denote the entropy measure in this 
joint distribution by M nma (P,Q). Then 

M nmia (P,Q) = M n , a (P) + M m , a (Q) + (a- l)M n>a (P)M m>a (Q). 

The third term on the right makes the measure non-additive. But 
the measures M% a (P) and Renyi entropy R k>a (P) as well as Shannon 
entropy S k are all additive. This is due to the fact that when logarithm 
of a product is taken it leads to a sum. 

(ix) Decomposibility: Consider a joint discrete distribution 

Pij > 0, Y^i=iYIj=iPij — 1- Consider the marginal distribution Pj = 
Y%=i Pij > 0, j — 1, . . . , m. Then we have 

M nm , a (pu,p 12 , . . . ,Pl m ,P21, • • • ,P2m, • • • ,Pnl, • • • , Pnm) 

= M m , a {p^ . . . , p m ) + jr pf- a M n , a ^). 

Observe that when p^ = piqj with Y%=iPi — 1> Y^JLiQj = 1 then 
property (ix) reduces to property (viii). 

(x) Functional equation: Consider M 2>a (P) = M 2;a (p, 1 ~ p)- That is, 

2-a 1 Q _ „)2-a _ 1 

M 2iQ (p,l-p) = ? ^ ^ , a^l, a<2. 

ct — 1 
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Let fa{p) = M 2t aip, 1 — p)- Then f a {p) satisfies the functional equation 

f a (x) + (1 - xf- a f a (^-) = f a (y) + (1 - yf- a f a (j^j (11) 

for x, y G [0, l),x + y G [0, 1] where [•) means open on the right and 
closed on the left, with 

/«(0) = f a (l) = (12) 

and 

/IN 2 a ~ 1 — 1 

Ms) = inm L < 13 > 

3.2 Characterization 

Several characterization theorems can be established for the measure of gen- 
eralized entropy in (5). We will indicate one such characterization. Let us 
look at an arbitrary continuous function, call it f a (x), satisfying the func- 
tional equation in (11) with the boundary conditions in (12) and (13). What 
is the functional form of f a (x)7 We will solve for f a (x) from (11). Put 
y = 1 — x in (11) and use the boundary condition (12) to obtain 

fa(x) = f a (l - X). 

Take two numbers p and q in the open interval (0, 1). Put p = 1 — x and 
q — y/(l — x) in (11). Then we obtain 

fa(p)+p 2 - a f a (q) = f a (pq) + (1 - pq) 2 - a fa 

Let 

F(p, q) = f a (p) + [p 2 ~ a + (1 - p) 2 - a ]f a (q). 
Then it is easily seen that F(p,q) = F(q,p), which implies that 

Up) + \p 2 - + (i - p) 2 -]Uq) = /«(?) + k 2 - + (i - q) 2 -°\U(p)- 



(IzJL) 

\l-pq) 
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Put q — 1/2 in this equation and use the boundary condition (13) to see that 



p 2-« + (l-pf-« 

fa{P) = 



a 



Now we can build up successively by using the branching property (viii). 
Theorem 2.1. Let fk(pi, ■■■,Pk) be an arbitrary function of pi, . . . ,pk, p% > 
0, i — 1, . . . , k, pi + . . .+Pk = 1 satisfying the properties of symmetry (prop- 
erty (Hi) of section 3.1), continuity (property (iv) of section 3.1), branching 
principle (property (vii) of section 3.1) where /2G0, 1 — P) — fa{p), satisfies 
the functional equation (11) with the boundary conditions in (12) and (13) 
then fkipi, ■ ■ -Pk) is uniquely determined as M k ^ a (P) of (5). 



3.3 Continuous analogue 

The continuous analogue to the measure in (5) is the following: 

MM) = sjumr^i (14) 

a — 1 

= : = : , « 7^1, «<2 

a — 1 a — 1 

where £"[•] denotes the expected value of [•]. Note that when a — 1, ^[/(rr)] 1- ' 
/•-• /(•'•) ° 1- 

When a < and decreases then 1 — a > 1 and increases. The measure of 
uncertainty decreases in the discrete case when a < 0. Similarly when ct > 0, 
then 1 — ct < 1 and decreases. In the discrete case the measure of uncertainty 
increases. Hence we may call 1-a as the strength of information in the dis- 
tribution. Larger the value of 1 — a the larger the information content and 
smaller the uncertainlty and vice versa. 



3.4 Connection to Kerridge's "inaccuracy" measure 

A connection to Kerridge's measure of "inaccuracy" can also be explored 
here. Kerridge (1961) defined 

X(P,Q) = -X>m%, (15) 
i=i 



8 



P= (pi, ...,p n ), Q = (qi, ...,q n ), Pi > 0, qi > 0, i = 1, . . . , n, £?=iP* = 
l,Z)iLi?i = 1, as a measure of "inaccuracy" in assigning the multinomial 
distribution Q for the true distribution P. Here Q may be relative frequencies 
coming from an experiment or observations and P may be the true underlying 
distribution. A generalization of (15) can be the following: 

M a {P,Q)= U = M , a?l, (16) 

a — 1 

so that when a — > 1 then M a (P,Q) goes to K(P,Q) of (15). Observe that 
we can also look upon M a (P,Q) as an expected value with respect to the 
true multinomial distribution P and we may write it as 

M a (P,Q) = E[ql a] ~\ a^l, (17) 
a — 1 

g taking the values qi, ■ ■ ■ ,q n with the corresponding probabilities Pi, ■ ■ ■ ,p n - 
The continuous analogue of (17) is the following: 

M a {f 1 ,f 2 ) = : = : , tt^ 1 (18) 

a — 1 a — 1 

where the expectation is taken with respect to the density fi(x). Instead 
of [/2(^)] 1_a if we assign [/i(a;)] 1_a through a displacement or distortion or 
disturbance to fi(x) then a measure of "inaccuracy" in assigning [/(x)] 1_a for 
f(x) is given by M a (f) of (14). More "inaccuracy" means less "information 
content" and more "uncertainty" or more "entropy", and vice versa. Again 
1 — a corresponds to the strength of information content. 

The generalized entropy M a (f) is evaluated for some standard distribu- 
tions and given in the appendix at the end of this paper. 



3.5 Distributions with maximum generalized entropy 

Among all densities, which one will give a maximum value for M a (f) ? Con- 
sider all possible functions f(x) such that f(x) > for all x, f(x) = outside 
(a, b), a < b, /(a) is the same for all such f(x), f(b) is the same for all such /, 
fa f( x )dx = 1. Let f(x) be a continuous function of x possessing continuous 
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derivatives with respect to x. Then for using calculus of variation techniques 
consider 

U = [f(x)r"-\ f(x). (19) 

[ b [f(x)] 2 - a dx-l 

Note that for fixed a, a ^ 1, maximization of — , , a ^ 1, a < 2 is 

equivalent to maximizing J^[f(x)] 2 ^ a dx. If necessary, we may also take 

Mo(/H _ slm±, a/1 , a<2 

a — 1 a — 1 

since J a f(x)dx = 1. This will produce only a change in the Lagrangian 
multiplier A in U above. Hence without loss of generality the form of U is 
as given in (19). We are looking at all possible / for every given x and a. 
Hence the Euler equation becomes, 

^ = => (2- aH/M] 1 -- A = 

free of x, a < 2, a ^ 1. Thus f(x) in this case is a uniform density over 
[a,b]. 

Let us consider the situation where E[x s ] for some 5 is a fixed quantity for 
all such /. Then we have to maximize 

f a [f(x)] 2 -dx 1 



a — 1 a — 1 

subject to the conditions f(x)dx = 1 and x 5 f(x)dx is a given quantity. 
Consider 

f/ = [/( a :)] 2 - a -A 1 /( a ;) + A 2a : 5 /(^). 
Then the Euler equation is the following: 

811 

= (2-a)[/(a;)] 1 - o '-A 1 + A 2 a; 5 = 

•9/ 

=► = ^-[1 - y^l 

2 — a Ai 

=S> /(x) = Ci [1-c 2 x 5 ]t^ (20) 
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where c\ and c 2 are constants and c\ > 0, 1 — c 2 x s > since it is assumed 
that f(x) > for all x. When c 2 = (3(1 — a), f3 > 0, we have 

f(x)= Cl [l-f3(l-a)x 6 ]T^. (21) 

Then for 5 = 1 we have the power law 

I = < 22 > 

where C3 is a constant. The form in (21) for a < 1 remains as a special 
case of a generalized type-1 beta model; for a > 1 it is a special case of 
a generalized type-2 beta model and when a — > 1 it is a special case of a 
generalized gamma model when the range (a, 6) is such that a = 0, b = 00. 
For 5=1, (21) gives Tsallis statistics (Tsallis, 1988, 2004). 
Observe that the generalized entropy M a (f) of (14) gives rise to the power 
law with exponent a, readily, as seen from (22). Also notice that by selecting 
Ai and A2 in (20) we can obtain functions of the following forms also: 

(1 - /^V 1 and (1 + (3 2 x s f\ 0!, fo,li,l2 > 0. 

Both these forms are ever increasing and cannot produce densities in (0, 00) 
unless the range of x with nonzero f(x) is finite. 

In section 3.4 we have given several interpretations for 1 — a. We can also 
derive the pathway model by maximizing M a (f) over all non-negative inte- 
grable functions. Consider all possible f(x) such that f(x) > for all x, 
fa f(x)dx < 00, f(x) is zero outside (a, b), /(a) is the same for all f(x), and 
similarly f(b) is also the same for all such functional /. Let f(x) be a con- 
tinuous function of x with continuous derivatives in (a,b). Let us maximize 
fa [f(x)} 2 ~ a dx for fixed a and over all functional /, under the conditions that 
the following two moment-like expressions be fixed quantities: 

f x^- 1)( - 1 - a) f(x)dx = given, and f 'a;(7-i)(i-«)+* f( x )dx = given (23) 

J a J a 

for fixed 7 > and 5 > 0. Consider 

U = [f{x)f- a - X lX ^- 1){1 - a) f(x) + X 2 x^- 1){1 - a)+s f(x), a<2, a^l 
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where Ai and A2 are Lagrangian multipliers. Then the Euler equation is the 
following: 

^ = (2-a)[/( a ;)] 1 - Q -A 1 ^- 1 )( 1 - a ) + A 2 ^- 1 )( 1 ^ +5 =0 
of 

=► [/(x)] 1 - = -^x^^ 1 -) [1 - ^V] (24) 
(2 — a) Ai 

=>• /(a;) = c 1 ^ 1 [l-/3(l-a)a; 5 ]T^ (25) 

where A1/A2 is written as — a) with /3 > such that 1 — (3(1 — a)x 5 > 
since f(x) is assumed to be non- negative. By using the conditions (23) we 
can determine c\ and (3. When the range of x for which f(x) is nonzero is 
(0, 00) and when c\ is a normalizing constant then (25) is the pathway model 
of Mathai (2005) in the scalar case where a is the pathway parameter. When 
7 = 1,5 = 1 then (25) produces the power law. The form in (24) for various 
values of Ai and A2 can produce all the four forms 

aix 7_1 [l - (3 X (1 - a)x 5 }~^, a 2 z 7_1 [l - (3 2 (1 - a)x s ]^ for a < 1 

and 

asx^ll + (3 3 (a - l)x 5 }-^, 04a: 7 " 1 [1 + (3 4 (a - l)x 5 }^ for a > 1 

with a i: Pi > 0, % — 1, 2, 3, 4. But out of these, the second and the third forms 
can produce densities in (0, 00). The first and fourth will not be converging. 
When f(x) is a density in (25) what is the normalizing constant c{l We 
need to consider three cases of a < l,a > 1 and a — > 1. This c\ is already 
evaluated in section 2. 

4 Reaction equation, superstatistics, and frac- 
tional calculus 

4.1 Reaction equation and fractional calculus 

A reaction problem was examined by Haubold and Mathai (1995, 2000) where 
the number density of the reacting particles is a function of time t. For the i- 
th particle let the number density be denoted by Ni(t) with Ni(t = 0) = Nq\ 
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If the production rate is proportional to the number density then we have 

j t Ni{t) = aiNi(t) 

where a. is the reaction coefficient. The reaction coefficient itself can be con- 
sidered to be a statistical quantity subject to accommodating a distribution 
function based either on Boltzmann-Gibbs or Tsallis statistics (Haubold and 
Mathai, 1998; Coraddu at al., 1999; Saxena, Mathai, and Haubold 2004b). 
If some of the particles produced are also destroyed and if the destruction 
rate is for the i-th particle, and proportional to the number density, then 

j t N % {t) = hN.it). 
If destruction dominates then the net residual effect is given by 

j t Ni{t) = -CiNi(t), c i = b i -a i >0 

and the solution is 

Ni(t) = N^e~ Ci \ a > o,f > o. 

But if a fractional integral, instead of the integer-order integral, is used then 
the residual reaction equation is given by 

Ni(t) = - c\ Dt v Ni(t), v > (26) 

where the fractional integral operator is given by 

a Drf(t) = ^ J\t - uy- l f(u)du, v > 0, (27) 

with a Dtf(t) = f(t). Note that in (26), q is replaced by c\ for physical 
reason. If the production rate dominates over the destruction rate then in 
(26), is to be replaced by — c\. For a discussion of general growth-decay 
models see Mathai (1993b). The Laplace transform of N^t) coming from 
(26) is 

s I 1 + \V . 
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and the solution can be written in terms of a Mittag-Lemer function: 

m = p = N ® Ev{ ~ <n (29) 

where the generalized Mittag-Lemer function E 1 (z) is defined by 
where (7) fe = 7(7 + 1)...(7 + A; - 1), (7)0 = 1, 7 ^ and 

00 k 



Consider the situation of having a total of 7 different particles, 7 = 1,2,..., 
these particles reacting independently with the same rate c\ — ci — ■ ■ ■ — 
c 7 = c. For independent variables the Laplace transform of the sum is the 
product of Laplace transforms. Hence in this case N(t) = Ni(t) + ... + N y (t) 
has the Laplace transform 

T («)- N ^--- N ^ (V) 
L N(t){s) - 77x1777 (olj 



1 + dr 



ST 

which yields the solution 

N(t) = . . . iV(S 7) F- 1 E%„(-f&'). (32) 

The rate c may vary and c may be considered to have its own distribution. 
In this case N(t) in (32) is to be taken as N(t\c) or N(t) at a given c. 
Superstatistics, as considered by Beck and Cohen (2003) and Beck (2006), 
will result when we use an appropriate prior density for c as developed in the 
following section. 

4.2 Superstatistics and fractional calculus 

Beck and Cohen developed the concept of superstatistics, see Beck and Co- 
hen (2003) and Beck (2006). They showed that superstatistics considerations 
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are relevant in problems of Lagrangian and Eulerian turbulence, defect tur- 
bulence, atmospheric turbulence, cosmic ray statistics, solar flares, and solar 
wind statistics. From a statistical point of view, the procedure is equiva- 
lent to starting with a conditional distribution for every given value of a 
parameter f3. Then j3 is assumed to have a prior known density. Then the 
unconditional distribution is obtained by integrating out over the density of 
P. 

Reaction equations are considered in a series of papers by Haubold, 
Mathai, and Saxena (Haubold and Mathai, 1995, 2000; Saxena, Mathai, 
and Haubold, 2004a; Mathai, Saxena, and Haubold, 2005). For the sake of 
illustration we consider one such reaction equation which yields the num- 
ber density of the following form, where 7 and /1 are arbitray, need not be 
integers. [Note that in (32) the resulting 7 is a positive integer.] 

N(t\c) = N Q t*- 1 Eff(-<?t v ), /i>0, 7>0, z/>0 

- N V (7+ !)*(-!) VP* m) 
~ Not to kWQt + vk) ■ (33) 

Let us consider the situation where c in equation (26) is a random variable 
having a gamma type density: 

g (c) = ^f-c^e-^, lo>0, 0<c<oo (34) 
1 (/ij 

where /x > 0,u,ii are known and /i/cu is the mean value of c. The residual 
rate of change may have small probabilities of it being too large or too small 
and the maximum probability may be for a medium range of values for the 
residual rate of change c. This is a very reasonable asumption. Equation (33) 
is the situation where the residual rate of change is such that the production 
rate dominates so that we have the form — c u , v > 0, c > 0. If the destruc- 
tion rate dominates then the constant will be of the form c u ,c > 0,u > 0. 
Integrating out N(t\c) over g(c) we have 



/ c =o ^ ' u r(/i) ^ h\r(ji + vk) 

where the integral over c yields r(/x + kv)oj~^ p+ku \ Hence the unconditional 
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number density, denoted by N(t), is given by 



= ^L t ^ 1 (l + -) for |-|<1. 

r(/i) V uo v ) V 

From the analytic continuation we see that the form remains the same for 
- > 1 also. Hence 

Nn / t v \ 

N{t) = Y^f~ l { 1 + ^J , < t < oo, cu> 0. (35) 

The continuation part may be seen by writing equation (33) with the help 
of the Mellin-Barnes representation in the form 

i f b+ioQ y(i — sirfs) , 

N(t\c) = N t^ — / I S) \ S) (c»tT s ds, i = V=T. (36) 

2lXl Jb-ioo 1 [l — sv) 

Then integrating over g(c) of (34) we have 

1 rb+ioo / +v \ s 

N(t) = N tr 1 — r(i- s )r( s ) - ds (37) 

2711 J6-ioo \UJ J 

where b is a real number, < b < 1. Evaluate the integral in (37) as the 
sum of the residues at the poles of T(s), which are at s = —k, k = 0, 1, to 
obtain the series form for £ < 1. Evaluate at the poles of T(l — s) to obtain 
the analytic continuation for £ > 1. Both will lead to the same form in (35), 
which is a generalized type-2 beta form. 

We may make the substitution 7 + 1 = — a > 1 =>- 7 = an d UJ ~ U — 
b(a — 1), b > 0. Then we have the unconditional number density 

N(t) = 7^ M-1 [1 + *(« " W 7 ^ (38) 

for a > 1, t > 0, b > 0, > 0. For /i — 1, v — 1, 6 = 1, equation (38) cor- 
responds to Tsallis statistics for one case, namely, a = q > 1 (Tsallis 2004). 
For general values of /i and a > 1 such that — ^ — ^ > equation (38) corre- 
sponds to the pathway model of Mathai (2005) as well as the superstatistics 
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considered by Beck and Cohen (2003) and Beck (2006). 

Now, consider an equation parallel to equation (17) in Mathai, Saxena, and 

Haubold (2005), namely, 

N(t\c) - Not^E-^n = c» D^N(t) (39) 

for 7 > 0, c > 0, v > 0, /j, > 0. Note that sets of parameters, 7 > 0, v > 
0,/i > exist so that the series form in E~^(c u t u ) is convergent and remains 
positive for all t > 0. This corresponds to the situation where the residual 
rate of change is positive so that the production dominates over destruction. 
We can expect N(t\c) to increase to 00. But the density of c in (34) will 
produce a dampening effect. Proceeding as before, we will obtain the final 
unconditinal N(t) as follows: 



for < - < 1. (40) 



For convenience, we can write 7 — 1 = -r—, a < 1 =>- 7 = f— - and uj v = 

1 1 l—a' ' I— a 

6(1 — a), b > 0, a < 1 so that (40) becomes 

N(t) = ^ )t ^[l-b(l-a)n^, (41) 

for a < 1,6 > 0,i/ > 0,pi > 0,0 < t" < [6(1 -a)]" 1 . If N(t) is to be made into 
a statistical density then we may normalize (38) with the help of a type-2 
beta integral, and (41) with the help of a type-1 beta integral. When a — > 1 
both (38) and (41) will approach the generalized gamma form 

N(t) = — ^"V 6 *" , 6 > 0, n > 0, v > 0, t > 0. (42) 

Equations (38), (41) and (42) are special scalar cases of the pathway model 
of Mathai (2005). For /i = 1, v — 1,6 = 1 Tsallis statistics is recovered from 
(38) and (41), with appropriate normalizations. 

Remark 1. Superstatistics considerations can produce only the generalized 
type-2 beta form of (35) with the appropriate normalization constant be- 
cause Beck and Cohen (2003) and Beck (2006) consider the situation of a 
conditional density belonging to the generalized gamma family where a scal- 
ing parameter is having a marginal distribution belonging to a generalized 
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gamma family. Then the unconditional density can only produce a form of 
the type 

f(x) = Cl x a (l + f3 1 x s )~ 1 \ ci > 0,/3i > 0, 7l > 0,x > (43) 

which is a generalized type-2 beta form. A type-1 beta form as in (40) will 
not be available from such a procedure because the exponents will only sum 
up to give a form of the type in (43). A more general format of the Beck- 
Cohen type consideration, from a statistical point of view, is provided in the 
following section. 

4.3 A more general density giving rise to generalized 
superstatistics 

Consider a conditional density of the form 

f x (x\(3) = c 2 (x - x ) a - l e- { ^ v{ - x ~^\ (44) 

for i] > 0, 5 > 0, a > 0, (3 > j3 Q , x > x Q where the normalizing constant c 2 can 
be seen to be the following: 

c 2 = . (45) 

r (f) 

For various values of the parameters a, rj, 5, x Q , (3 we have the exponential, 
gamma, chisquare, Erlang, Helley, Helmert, Maxwell-Boltzmann, Rayleigh, 
and Weibull densities as special cases of (44) (Mathai, 1993a; Honerkamp, 
1994). If x — x is replaced by \x — x \ for — oo < x < oo then many more 
densities such as double exponential, Laplace and Gaussian will be particular 
cases. Consider a prior density for the scaling parameter f3 in the following 
form: 

g{P)=<%{p-p o y- 1 e<W-M\ (46) 

for ( > 0, rj > 0, 7 > known, (3 > (3 and c 3 = /T (j^. This is a very 
general class of the type in (44). The only restriction imposed for convenience 
is that the exponents involving (3 are the same in (44) and (46). Then the 
unconditional density of x is given by 



fx(x) = _ f x (x\(3)g((3)d(3 

Jj3—(3 
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(/3-/3 o ri(+(x-Xo) s ] d p 



C f r(j)r(f) 




(47) 



for x > x , ( > 0, 7 > 0, i] > 0, a > 0, 5 > 0. Note that (47) can only produce 
a form of the type (43) or a generalized type-2 beta type. Type-1 beta form 
of (40) is not available from a procedure such as the steps in (44) to (47). 
This is not covered in Beck-Cohen's type superstatistics considerations. 
Remark 2. Observe that if a real random variable x has a density belonging 
to the generalized gamma family then ^ also belongs to a generalized gamma 
family. Hence in (44) and (46) we could have taken (x — x )~ 5 with 5 > and 
(f3 — Po)~ v with rj > in the exponents, individually or simultaneously. The 
procedure will remain the same and we will again end up with a generalized 
type-2 beta form corresponding to (47). 
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Appendix 



Al. Evaluation of the generalized entropy for some standard dis- 
tributions 



Gamma density: 



for x > 0, 7 > 0, (3 > and f(x) 
and exponential densities also. 



elsewhere. This includes the chisquare 



/•oo 

/ [f(x)] 2 - a dx = 
Jo 



Jo x 



( 7 -l)(2-a), 



(2-a)j: 



/3 dx 



/3(2-ah[r( 7 )] 



2-a 



r- 1 r[( 7 -l)(2-a) + l] 



[r(7)l 



2-a 



Hence 



M Q (/) 

and 
limM a (/) 

a— >1 



r- 1 r[( 7 -l)(2-a) + l] 

" [r( 7 )] 2 - 



ln/3 + (l -7)^(7) -lnr( 7 ) 
In /3 for 7 = a, 



where ip(-) is a psi function or digamma function. It may be also noted that 
the continuous analogue of M% (P) of (9) is 



K(f) 



a 



1 



-, a 7^ 1, a < 2. 



For the gamma density M*(f) also gives the same quantity as above, and 
In (3 when a — > 1 . 

It may be observed that the same procedure goes through even for the path- 
way model of Mathai (2005), in the scalar case, and M a (f) as well as M*(f) 
can be computed explicitly without much difficulty The pathway model con- 
tains a large number of distributions of common use in various disciplines. 
But Shannon entropy Sk of (6) is quite difficult to evaluate, in the continuous 
analogues, even for many of the standard densities in common use due to the 
presence of In f(x). 
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